Decohering the Fermi hquid: A dual approach to the Mott Transition 
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We present a theoretical approach to describing the Mott transition of electrons on a two dimen- 
sional lattice that begins with the low energy effective theory of the Fermi liquid. The approach to 
the Mott transition must be characterized by the suppression of density and current fluctuations 
which correspond to speciflc shape deformations of the Fermi surface. We explore the nature of the 
Mott insulator and the corresponding Mott transition when these shape fluctuations of the Fermi 
surface are suppressed without making any a prior assumptions about other Fermi surface shape 
fluctuations. Building on insights from the theory of the Mott transition of bosons, we implement 
this suppression by identifying and condensing vortex degrees of freedom in the phase of the low 
energy electron operator. We show that the resulting Mott insulator is a quantum spin liquid with a 
spinon fermi surface coupled to a f/(l) gauge fleld which is usually described within a slave particle 
formulation. Our approach thus provides a coarse-grained treatment of the Mott transition and 
the proximate spin liquid that is nevertheless equivalent to the standard slave particle construction. 
This alternate point of view provides further insight into the novel physics of the Mott transition 
and the spin liquid state that is potentially useful. We describe a generalization that suppresses spin 
anti-symmetric fluctuations of the Fermi surface that leads to a spin-gapped charge metal previously 
also discussed in terms of slave particle constructions. 
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PACS numbers: 



I. INTRODUCTION 



Describing the evolution from a metal to a Mott in- 
sulator as the strength of electron correlations is in- 
creased has challenged condensed matter theorists for 
the last many decades. Part of the difficulty is due to a 
complicated interplay between the metal-insulator tran- 
sition and magnetic ordering that appears commonly in 
the insulator. However in recent years evidence, both 
theoretical— and experimentalSr-, has accumulated for 
the existence of a possible quantum spin liquid Mott state 
just on the insulating side of the Mott transition in frus- 
trated lattices. Such a state does not break spin rotation 
or any other microscopic symmetries. Thus there is an 
opportunity to study the fundamental phenomenon of 
the Mott metal-insulator transition without the compli- 
cations of magnetism or other broken symmetries. 

A theoretical description of the vicinity of the Mott 
transition is challenged by the need to describe the ki- 
netic and interaction energies of the electrons on more 
or less equal footing. In the last two decades there has 
been great progress in studying the Mott transition in 
spatial dimension d = oo through Dynamical Mean Field 
Theory^*'. However in this limit the resulting Mott insu- 
lator has at zero temperature no correlations between the 
spins and hence has non-zero entropy. Thus the d = oo 
results may at best be a guide for the somewhat high 
temperature physics in the vicinity of the Mott transi- 
tion. 

The Mott transition can also be described using slave 
particle methods. A complete theory for a continuous 
Mott transition from a Fermi liquid metal to a quantum 
spin liquid Mott insulator in two space dimensions was 
described in Ref. [ll| by studying fluctuations on top of 
the slave particle mean field theory of Ref. [l^- In this 



description the Mott insulator is a quantum spin liquid 
state with gapless spin excitations. As this transition 
is approached from the metal the entire Fermi surface 
disappears continuously. The resulting critical point is 
characterized by the presence of a gapless 'critical Fermi 
surface— where the Fermi surface is sharp but there are 
no Landau quasiparticles. These results have also been 
generalized to three space dimensionsi^. 




FIG. 1: A: Fermi surface at equilibrium. B: The low energy 
excitations describe smooth shape fluctuations of the Fermi 
surface. 

In this paper we develop a way of thinking about the 
electronic Mott transition in two space dimensions that 
starts with the low energy effective theory of the Fermi 
liquid phase. This is of course just the familiar Landau 
theory of Fermi liquidsi^. The hydrodynamic descrip- 
tion of the low energy properties of a Fermi liquid is in 
terms of the kinetic equation for the local quasiparticle 
distribution function. This hydrodynamic theory may 
be usefully viewed as a theory of shape fluctuations of 
the Fermi surface^^"— (see Fig. [Ij. In the low energy 
limit a local version of Luttinger's theorem22 is obeyed 
so that a local shape fluctuation of the Fermi surface 




FIG. 2: A: Current fluctuation {I = 1) of the Fermi surface. 
B: Density fluctuation {I = 0) of the Fermi surface. 



which changes the area leads to a change of particle den- 
sity (see Fig. [2]). Similarly a local shift of the "center 
of mass" of the Fermi surface corresponds to a non-zero 
current density (see Fig. [5]). As the Mott transition is 
approached the minimal thing that needs to happen is 
that long wavelength fluctuations in the local charge and 
current densities are suppressed. Our goal in this pa- 
per is to explore the nature of the Mott insulator and 
the transition when the Fermi surface shape fluctuations 
associated with charge and current density fluctuations 
are suppressed without making any a priori assumptions 
about the fate of other shape fluctuations. 

We will do this by relying crucially on key insights from 
the theory of the Mott transition of bosons. In contrast 
to electronic systems, the Mott transition of interacting 
bosons is rather well understood theoretically^^. The 
liquid phase of bosons has superfluid long range order, 
and the corresponding hydrodynamic theory describes 
the gapless sound mode associated with superfluid order. 
Accessing a Mott insulator of bosons from the superfluid 
phase requires condensing vortex configurations of the 
superfluid order parameter. This enforces the quantiza- 
tion of particle number that is crucial for the existence 
of a Mott insulator. 

Motivated by this, we develop an alternate 'dual' way 
of thinking about the electronic Mott transition and the 
resulting spin liquid Mott insulator. We first identify 
a degree of freedom in the hydrodynamic description of 
the conducting state that plays the role of a vortex in the 
phase of the electron operator. We show that condensing 
these vortices leads to an incompressible Mott insulator. 
Specifically it clamps down long distance, long time fiuc- 
tuations of the particle number and current densities. 

The resulting Mott insulator is in a quantum spin 
liquid state and has gapless excitations. Interestingly 
it is a spin liquid with a Fermi surface of charge neu- 
tral spin-1/2 fermionic spinous coupled to a fiuctuating 
C/(l) gauge field. Usually spin liquids with emergent 
spinous and gauge fields are described by slave particle 
constructions^i. The route to a spin liquid described in 
this paper is very different and does not explicitly rely 
on any slave particle description. It thus provides an al- 
ternate point of view of the spinon Fermi surface state 
that may potentially be useful in thinking about it. We 



also briefly discuss a phase obtained by condensing vor- 
tices in a phase that is conjugate to spin antisymmetric 
distortions of the Fermi surface. This results in a spin 
gapped metal with a charge Fermi surface which can also 
be alternately described in terms of slave particles. 



II. OVERVIEW 

A. Bosons and the Mott transition 

Consider, for concreteness, the boson Hubbard model 
on, say, a two dimensional lattice with an integer number 
iVo of bosons per site on average: 
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where b] creates a boson at site i and Ui = b\bi is the bo- 
son number at site i. liU ^ t the ground state is a Mott 
insulator while in the opposite limit t ^ U a. superfluid 
state results. The superfluid-Mott phase transition is 
describec^'^ by a standard (quantum) Landau-Ginzburg 
action for the superfluid order parameter in Z? = 2 + 1 
space-time dimensions. 

For our purposes it is insightful to understand how 
to think about the Mott insulator and the transition to 
it within in a conceptual framework that begins in the 
superfluid phase. Deep in the superfluid state the appro- 
priate low energy effective theory of the superfluid phase 
is simply given by the harmonic (Euclidean) action 
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(2) 



for the phase cf) of the superfluid order parameter (6 



). Here ^ = {^tt) and we have set the sound velocity 



ai<t> 

to 1. ps is the phase stiffness. A key point is to recognize 
that this "phase-only" description has no hope of describ- 
ing the Mott insulator. This is because the harmonic 
phase-only theory does not know about the quantization 
of the conjugate particle number. This quantization is 
crucial in obtaining the Mott insulator. Quantization of 
the particle number is equivalent to the condition that 
the boson phase be defined mod 2-k. Thus it is impor- 
tant to include vortices in the phase field and condense 
them to obtain a description of the Mott insulator. As is 
well known'^^''^^ this is conveniently described through a 
duality transformation where we regard the gapless lin- 
early dispersing sound wave as a transverse photon of a 
C/(l) gauge field. Formally 



5cff -^ / (fxdr (V X a) 



(3) 



Vortices in (/) act as sources for the dual vector potential 
and can be easily incorporated as a boson field that 



^/i 1 



couples minimally to a^. The full dual action including 



a vortex field <J>i, is 

S^ f d^xdT^ (V y.af + \ {d^ - la^) $„|2 + V[\^^\^) 

(4) 
In the superfluid phase the vortices are gapped and may 
be integrated out to reproduce the dual of the low energy 
sound wave action. The Mott insulator is obtained when 
the vortices condense. This gaps out the gauge field, and 
quantizes the dual magnetic flux which is the particle 
number of the original bosons. The dual vortex action 
is well known to be equivalent to the Landau-Ginzburg 
description of the Mott transition in terms of the boson 
order parameter h. 



B. Electrons and the Mott transition 

So now consider the electronic Hubbard model at half- 
filling on a non-bipartite lattice in two space dimensions: 
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where c\^ {a =ti i) creates a spinful electron at site i 
and Tii = ^^ CmCia is the electron number at site i. For 
t ^ U a Fermi liquid results with a sharp Fermi surface 
satisfying Luttinger's theorem. For U ^ t on the other 
hand a Mott insulator obtains. At very large U the low 
energy physics of the Mott insulator may be described 
as a spin model and will often have magnetic long range 
order. However we will specifically be interested in the 
'weak' Mott insulator that obtains just on the insulating 
side of the Mott transition. It seems likely, for instance 
on a triangular lattice, that such a weak Mott insulator 
will be in a quantum spin liquid state with no magnetism 
or other broken symmetriesii^ii^iii^. 

It is interesting, in the fermionic version of the prob- 
lem, to start again with the low energy effective theory of 
the conducting phase and understand what needs to be 
done to access a Mott insulator. In the fermionic case, 
as the conducting phase (for t ^ U ) is just the Fermi 
liquid, the corresponding low energy effective theory is 
Fermi liquid theory. Analogously to the bosonic example, 
the Landau Fermi liquid action by itself does not contain 
within it the Mott insulator that obtains with increasing 
U. We may then ask if there are special configurations 
(analogous to the vortices in the boson example) that 
have to be included and whose condensation can lead to 
the Mott insulator. The analogy with bosons may be 
sharpened somewhat by noticing that Fermi liquid the- 
ory is essentially a quadratic hydrodynamic theory for 
the density of fermions at each point of the Fermi sur- 
face. (For bosons the hydrodynamic mode is just the 
conserved total density and the quadratic phase theory 
just describes fluctuations of this mode. In contrast for 
fermions the Fermi liquid 'fixed point' has emergent con- 
servation laws corresponding to the number of fermions 



at each point of the Fermi surface^^, and the correspond- 
ing densities are emergent hydrodynamic modes that are 
included in the Landau Fermi liquid description). Fur- 
thermore this is equivalents*'^-' to a quadratic 'phase' 
based description in terms of the phase of the fermion 
operator at any given point 9 of the two dimensional 
fermi surface. This quadratic phase action is a higher 
dimensional generalization^^^^^ of the familiar bosoniza- 
tion from one spatial dimension. However in contrast to 
the one dimensional example the bosonized Fermi sur- 
face theory in higher dimensions is simply a rewrite of 
Landau fermi liquid theory. 

For the purposes of accessing a Mott insulator, it 
seems most interesting to explore a slightly different 
point of view on the Fermi liquid. As emphasized in 
the literaturei^iii, the quadratic hydrodynamic theory 
for the emergent conserved densities of the 0- movers (i.e. 
fermions associated with an angle on the Fermi surface) 
is a theory of harmonic fluctuations of the shape of the 
Fermi surface (see FIG.[T]). 

The fluctuations of the total density are simply deter- 
mined by the fluctuations in the area of the Fermi surface 
(with the assumption^^ that long wavelength fluctuations 
locally satisfy Luttinger's theorem). Clearly to access the 
Mott insulator a description where the quantization of 
particle number is incorporated is necessary. By analogy 
with bosons we may then expect that we need to include 
vortices in the phase variable that is conjugate to the 
total particle density. This conjugate phase variable is 
readily identified in the bosonized phase description of 
the Fermi liquid (as a suitable linear combination of the 
phase of the various 0- movers). We may then anticipate 
that a modified theory that correctly includes vortices 
in this phase variable will enable description of a Mott 
insulator. 

The rest of the paper is structured as follows: In Sec. 
mil we will briefiy review bosonization of a Fermi liquid 
in 2+1 dimensionsi^ii^. In Sec. IIVI we show how to in- 
corporate vortices into the phases of the fermions and 
thereby obtain an effective 'dual' low energy theory that 
generalizes Fermi liquid theory. In Sec. |V] we discuss 
the phase structure of this dual theory and show that it 
contains both Fermi liquid and Mott insulating phases. 
Furthermore we show that the Mott insulator is a quan- 
tum spin liquid that is described at low energies in terms 
of a spinon Fermi surface coupled to an emergent U{1) 
gauge field, in agreement with more standard slave par- 
ticle approaches. Then in Section I VII we show how by 
condensing a vortex in the phase conjugate to spin an- 
tisymmetric fiuctuations of the Fermi surface, we may 
access a non-fermi liquid metal with a spin gap but a 
charge fermi surface. In Section I VIII we briefly consider 
defining and condensing vortices in phases conjugate to 
shape fiuctuations of the Fermi surface with non-zero an- 
gular momentum. Section. IVIIII contains a summary of 
our results. 



III. BOSONIZATION OF A FERMI LIQUID IN 

2D 



Landau Fermi liquids possess a sharply defined Fermi 
surface where the electron momentum distribution has 
a jump discontinuity. The low energy excitations may 
be described in terms of an electron-like quasiparticle 
that is well defined asymptotically close to the Fermi 
surfacei^^. The long distance, low energy properties 
of the Fermi liquid may be described in terms of a quasi- 
particle distribution function n„(r,k,t) (with the spin 
index a =t, 4-) which is well defined for k close to the 
Fermi surface, and when 5na^k[x,t) = na^k{x,t) — ?^" ^ 
is small, n" = 0(/x — e^) is the ground state distribution 
(/U is the chemical potential and e^ is the dispersion of 
a single quasiparticle). A configuration 5n{r,k) has the 
energy 
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(6) 
Here V is the system size and /^X' i^ the familiar Lan- 
dau interaction function. In the low energy limit we can 
ignore quasiparticle collisions and the 5n satisfy a (lin- 
earized) kinetic equation 
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where vp = ^e^ and terms of 0{Sn^) were dropped. 

To address the low energy, long wavelength physics we 
restrict to a small band of momentum width A|| near the 
Fermi surface, and define 
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Here 9 is an angle that labels points on the Fermi sur- 
face and the corresponding Fermi momentum is Kp{9). 
The integration variable k is the deviation of the mo- 
mentum in the direction normal to the Fermi surface at 
point 6. Clearly we may interpret A6Kp{9)pcr{9,r,t) as 
the density of electrons (of spin cr) in a patch of angu- 
lar width A9 centered at point 9 of the Fermi surface. 
We will dub these electrons 6'-movers by analogy with 
the familiar terms left /right movers in d = 1 Luttinger 
liquids'^^. The deviation of the total electron density 
6p^(r,t) from its mean is clearly related to the density 
of 0-movers through 



(9) 
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It will sometimes be useful to work with discrete 
patches: chop up the Fermi surface into N patches of an- 
gular width A9 
we may write 



The 0-mover density may be given an interpretation in 
terms of the fluctuation of the shape of the Fermi surface. 
Consider a slow long wavelength distortion of the Fermi 
surface where 



Kpa ^ Kpa + SKpair, 9, t) 



(11) 



Assume that the Luttinger theorem relating the area of 
the Fermi surface to the density holds locally even for 
such a long wavelength disturbance of the Fermi fluid. 
The the change in the density is 



5 Pair, t) 
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The change in the density of 0-movers may then be iden- 
tified as 



p{r,9,t)^^6KFa{r,9,t) 



(13) 



We may now reduce the full kinetic equation to write 
down the equation of motion of the 0- movers. To that end 
we first let vp,the the normal and tangential unit vectors 
associated with point 9 on the Fermi surface. Then by 
integrating with respect to fcy we find 

dtP.ie)+VFd^^p4e)+j^ j de'KF{e')u^,{e, e')p,,{e') = o 

(14) 
Here d\\ = vf ■ V is the derivative in the direction normal 
to the Fermi surface. 

A bosonized phase representation is obtained by intro- 
ducing for each patch i a phase field 0^0- such that 



Pia = -z-^\\4>i<y 
111 



(15) 



Substituting in Eqn. [T3] we get an equation of motion for 
(pia which takes the form 

(9t +Wf5||) ^\\(j,^a + A9Y,KF^' Ua'{hi')d\\d\\(j)^,a' = 

(16) 

This is clearly a direct generalization of the familiar 
phase representation of the right and left moving densi- 
ties in one dimensional liquids. Just as in one dimen- 
sion (pia will be interpreted as the phase of the ^-moving 
fermion, as we now discuss. The equation of motion 
above can be obtained from the Lagrangian 

L = Lo + Lf (17) 

A9 



;ie reiuii Buriace niLo jv paicnea ui an- AW -r-^ / / x2\ 

2f and eventually let N ^ oo. Then ^o = —^^^Ft [-^t(|}^a^\l(|3^a - vf (5|| </),,,) j (18) 
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Lf = - 
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We can now write the quantum partition function as an 
imaginary time path integral 



(20) 
(-idT(l>i^d\\(j>ia+VF{d\\(j)ia) ] (21) 
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The structure of the time derivative term in the La- 
grangian determines - in an operator framework - the 
commutation relations for (f>i^. Of crucial importance is 
the following relation shown in the Appendix \X\ 
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Here {x,y)\\ are the components of the spatial coordi- 
nates in the direction parallel to the patch normal and 
{x, y)± are the components in the direction perpendicular 
to the patch normal. The function Z? is a delta-function 
smeared over a distance of order the inverse patch size 
^ K Ae ■ Thus in the long wavelength limit 2TT(j)i is con- 
jugate to the density of ^-movers. We may now sum over 
all patches to find 
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with (j)Qcr defined by 
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Thus 27r0O(T is canonically conjugate to the long wave- 
length fluctuations of the total density of spin a. Defining 
the angular Fourier transform of the <j)i fields: 
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we identify ^oo- with the / — component. It is thus 
natural that it is conjugate to the uniform "breathing" 
mode of the Fermi surface which corresponds to a change 
of electron density. If one wishes the quadratic action for 
the (t>ic fields can be readily rewritten in terms of Fourier 
transformed variables (pia- 

The commutation relation Eqn. [531 implies that the 
operator e^^"^'" adds charge 1 to patch i with spin a. 
In the Appendix VK\ we further show that we may indeed 
identify this with the electron creation operator at patch 
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As in one dimensions it is now convenient to go to 
charge and spin bosons (j)ci,4'si- We thus define 



2 
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The action now splits into separate ones for (pd, 4>si'- 

S = ScQ + Scf + SsQ + Ssf (30) 

ScO=J2 I ^^^^» (-i^r(|>c^^\\<|}o^ + Vp (9||<?icO') (31) 
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Here we defined the standard symmetric and antisym- 
metric combination /g, fa of the Landau interaction func- 
tion (/, = ^^^^^Ja = ^^^^^). The commutation re- 



lations satisfied by 
Clearly we have 



may be readily written down. 
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{x),Sp{y)] = ^S^'\x-y) (35) 



with the obvious notation i^co = 2 °^ ' ^P — '^a ^P^ 



IV. INCORPORATING CHARGE VORTICES 

The Mott insulator is incompressible for density fluctu- 
ations. Thus we need to clamp down fluctuations in the 
shape of the Fermi surface that correspond to a change 
of area. Furthermore we also need to clamp down fluc- 
tuations where the center-of-mass of the Fermi surface is 
displaced - these correspond to a non-zero current density 
(see Fig. H]). 

Following the logic for bosons we expect that these 
can be accomplished by including vortices in the phase 
of the fermions. Specifically if we include vortices in the 
phase (/iQc (which is conjugate to the total density) and 
condense them, we will get a fermionic Mott insulator. 

To think about a vortex state in a Fermi liquid, first 
consider a ring geometry obtained by imposing periodic 
boundary conditions along the a;-direction, i.e. we iden- 
tify X with X + L. Imagine slowly turning on a 27r flux. 
This has the effect of twisting the electron boundary con- 
ditions on going around the ring by 2tt. Correspondingly 
all single particle momentum eigenstates shift by ^^. 
For a free fermi gas the ground state wavefunction is 
then given by 

i^twi.t)= n 4+2j.A^,iio) (36) 

\k„\eFS 

This shifts the entire Fermi surface by ^^. Therefore 
the Fermi surface displacement is 



2n 
SKpaid) = -j-vpe ■ X 
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which can be obtained from 
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(37) 
(38) 



This is consistent with the expectation that the phase 
of the electron is shifted by 2tt when a 2tt twist of the 
boundary condition is imposed. Clearly 



dl ■ V<j)^a = 1 



(39) 



for any loop around the ring. Thus there is a 27r vortex 
in the electron phase. This vortex is present in both (p^ 
and (j)^. Thus (pc winds but (f>g does not wind at all. 

We may also directly consider such vortices in the bulk. 
Consider a "core" such that the electron phase (pco winds 
by 27r on going around a loop surrounding the core. In 
the presence of such a vortex (pc will not be a smooth 
field anymore. We now show how the Fermi liquid ac- 
tion may be generalized to allow for the presence of such 
vortices. As in the usual discussion of the 2d XY model, 
we separate 



•co into a smooth part (p and a vortex part 
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The smooth part satisfies 



e^-A^.^A^^o = (41) 

The vortex part is defined in terms of the vortex 3-current 
j'^i, through the equation 



f^t^u\dud\cpl^^ = ji_,, 



(42) 



As 5^(/)co is not really a simple gradient we will replace 
it by a vector field a^ satisfying 



^^u\^u(^\ — 3 ^lv 
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Thus vortex degrees of freedom may be included in the 
Fermi liquid action by simply replacing every occurrence 
of df^(pco with df^(p „ + a^. Some care is required however 
in implementing this procedure. To see this consider the 
terms in the bosonized action involving time derivatives: 
it has the structure 



J T.r 



(44) 



with Ki = KpiVpi- In terms of the Fourier modes (pd, 
this becomes 



/ y^ -~idr(pclK_i_i,.V(pcl' 



(45) 



with Ki = ;^ X^i ^~''^'-^i- -fo'" smooth configurations 
of all the (pi fields, this action is symmetric under inter- 
change of I and /' - we simply integrate by parts twice to 
also interchange the time and spatial derivatives. This 
can be exploited to rewrite the action in several equiva- 
lent forms - for instance we may symmetrize with respect 
to / and V . However for non-smooth configurations of cpi 
these different forms are not equivalent, and we need a 



prescription to choose between them. The physically sen- 
sible prescription which we will employ is to insist that 
(pQc enters the time derivative term in the action in such 
a way as to ensure that it correctly obeys- in an operator 
framework- its commutation relation with the total den- 
sity. To implement this, we separate terms that involve 
(poc from the rest to rewrite Eqn. |45] as 



/ -2idr<PcoYK^vV(pcl'-\ ^ dr<PclK^i^i>-V(pcl' 

(46) 
In the first term the sum over I' does not receive any 
contribution from /' = as 
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This follows from time reversal invariance or inversion 
symmetry. Then (pco enters this part of the action linearly 
and only through its time derivative. Therefore we can 
directly read off its canonical conjugate and check that 
it is indeed the total density. We find 
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which is readily seen to agree with our earlier expression 
for 6p. 

As we are concerned with vortices in 0co and none in 
any other (pci,l ^ 0, why do we not simply work with 
the (pci for each patch and separate out a common vortex 
part? The answer has to do with the ambiguity men- 
tioned above. The same issue arises when we try to cou- 
ple an external gauge field to the bosonized Fermi liquid 
action, and in the theory of d = 1 fermionic Luttinger liq- 
uids. Indeed in the presence of an external gauge field, 
separate conservation of the densities of 9 movers is vio- 
lated as the entire Fermi surface is displaced by the elec- 
tric field. This is the famous anomaly"^"'^ familiar from the 
theory of fermions in d — 1. This means that we are not 
free to minimally couple a gauge field or to extract the 
vortex part independently from each (pd ■ Rather the cor- 
rect procedure is indeed to work with 0co and proceed as 
we did above to either introduce vortices or to couple in a 
gauge field. We illustrate this issue further in Appendix 
|B]in the familiar context of d = 1 Luttinger liquids. 

We may now include vortices in the bosonized action 
by separating out the smooth part of d^<pco from its vor- 
tex part as discussed before. For notational convenience 
we shall drop the underscore from the smooth part from 
now on. The resulting Lagrangian density for the charge 



boson then takes the form 
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Here M^, is the angular Fourier transform of the kernel in 
the terms in the original charge action that depend only 
on spatial derivatives. We note here that the vortex part 
a^t couples in to the charge boson in exactly the same 
form as a U{1) gauge field A^'^^'' , which will be impor- 
tant to us below. To see this it is sufficient to consider 
^gauge ^ ^gaugo ^ Q, as the ambiguity discussed above 
(see also Appendix [B]) only involves the time-derivative 
term in the action. The temporal component of the gauge 
field Af^"sc couples to the total charge density, i.e. 



Cc{At 



augc 



Cc 



ipAf "g'' . (53) 



■'C I J^S^MgC — Q 

But this is identical to ^ for Af-'^SB _^ 2'Kar, where 
the 2tt reflects the identification of 27r(^9 as the phase of 
a fermionic 0-niover. 

Now we introduce a vortex field $„ whose three-current 
is precisely j^^ (see eqn. |35]). Then the modified charge 
Lagrangian density is 

£ = Cc[d^(j)ci + afj,5ifi] + iA^€fj,i,xd^ax + £[$„, A^] (54) 

where A^ is minimally coupled to $„. Varying with re- 
spect to Af^ then implements Eqn. l42l 

We now specialize to the situation that the total charge 
density of electrons is such that there is one electron per 
site on average. Then the average flux seen by the vor- 
tices will be zero, and the vortex Lagrangian will have 
the structure 



£[$,, Aj = I (a^ - iA^) $„|' + y (1$. 



) + 7^-2 (i^t^i^xdiyAx) 
ZCq 

(55) 
where y is a potential which can be expanded as a poly- 
nomial. 

The resulting full charge Lagrangian described by 
Eqns. |49l [54l and [55] then taken together with the spin 
Lagrangian completes our modification of the basic Fermi 
liquid action to incorporate charge vortices. 



V. PHASE STRUCTURE OF DUAL VORTEX 
THEORY 

We will now use this dual vortex theory to discuss var- 
ious possible phases of the fermion system. 



A. Fermi liquid 

When the vortices are gapped we expect to recover the 
Fermi liquid. To see this explicitly, we integrate out the 



$« field from the action. In the presence of the vortex 
gap, this will lead to an innocuous renormalization of the 
coefficient of the Maxwell term for the A^ field: 

1 2 



(56) 



Now integrate out A^ 



VA Q^L,-r'^^l'<^l"'X^'^°-^ + ^^'^t''^>-^"^ 



-/.. 



(57) 

Thus the a^ get a 'mass' term and can be ignored for 
the low energy physics. Setting a^ = we obtain just 
the bosonized action for the Landau Fermi liquid that 
we started out with. 



B. Mott insulator 

When the vortices condense the fluctuations of the 
fermion density and current will be suppressed and we 
will get a Mott insulator. Our formalism will enable 
us to obtain a 'dual' description of this Mott insulator. 
When ($„) ^ 0, there will be a Meissner effect for A^, 
and it will acquire a 'mass'. Integrating out Ap, will now 
generate a Maxwell term for a^. Then the full effective 
Lagrangian in the Mott insulating phase 

K 2 

£cff = Cc[df,(l)ci + a^(5;,o] + Y {e^j^uxduax) (58) 

where ii' is a non-universal coefficient. Rewriting this in 
terms of 4>^ , 4>i we get 



£cff = £[9p(/)jt + a^'^/,0, 9^0-it + a,^,6ifi] 



y {etj.uXOuax) 

(59) 
This is the same as the Fermi liquid action except that the 
phases 0^, (jy^ are minimally coupled to a C/(l) gauge field. 
This observation enables us to now go back to a fermionic 
representation. We introduce fermion fields fai ^ e'^^'^^" 
for each patch i and spin a. The structure of the low 
energy effective theory for the Mott insulator is simply 
that of fermions at a Fermi surface coupled minimally to 
the fluctuating gauge fleld a^. The corresponding action 
may be taken to be 



m, Am] = faidr - iao - f-i - 



Y (f^^LuxOvax) 



{V-iaY 
2m 



fa 



(60) 



Clearly the /„ field should be interpreted as a spinon, 
and the Mott insulator we obtain is described at low en- 
ergies as a spin liquid with a spinon fermi surface that is 
coupled minimally to a fluctuating U{1) gauge field. This 
is the same result as that obtained within the standard 
slave particle treatment of the Mott transition^i. How- 
ever we have reached it through a very different route 
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that never involved splitting the electron into a product 
of slave particle operators. Thus our approach provides 
an alternate point of view on the spinon fermi surface 
spin liquid Mott insulator that adds to our understand- 
ing of it. 

Let us pursue a bit more this connection between the 
approach we have taken in this paper and the standard 
slave particle description of this spin liquid Mott insula- 
tor. We see that the chargon/holon operator of the slave 
particle formalism should be identified with e'^^^'^"", i.e. 
the holon/chargon phase is essentially the I — compo- 
nent of the charge boson. The bosonized version of the 
spinon field with spin a at angle 6 on the Fermi surface 
is simply g2i"(0<T(e)-</>oc)^ Finally the magnetic flux of the 
emergent gauge field of the spin liquid is simply the den- 
sity of the vortex field $„ . 

In the slave particle approach the internal gauge field 
is compact. The compactness means that "instanton" 
events where the gauge flux changes by 2tt are microscop- 
ically allowed. It has now been established that in two 
dimensional spin liquids with sufficient number of gapless 
spinon modes these instanton events are supressed at low 
energies so that the compactness can be ignored in the 
low energy effective theory of the spin liquid phase^^^. 
How does this manifest itself in our approach? Since the 
magnetic flux of the emergent gauge field corresponds to 
the density of vortices, these instanton events correspond 
precisely to processes where the vortex decays and dis- 
appears. This is exactly what one expects in the Fermi 
liquid phase. As we discussed a vortex corresponds to a 
pattern of swirling current and in the Fermi liquid phase 
it will decay into a shower of particle-hole pairs. This 
decay process of the vortex is actually not included in 
our effective dual action. However as we know that the 
instanton events are suppressed in the spin liquid phase 
and right at the Mott transition we are justified a poste- 
riori in ignoring this vortex decay. 



VI. SPIN- VORTEX CONDENSATION AS A 
ROUTE TO A SPIN GAPPED METAL 

In this section we briefly digress and consider the nat- 
ural question of what happens if we condense the vortex 
in the spin rather than in the charge boson. Specifically 
let us define the zero angular momentum component of 
the spin boson: 



«o-]^E 



(61) 



It is easy to see that tt^cO is conjugate to (twice) the 
^-component of the spin density M^ — (Sp^ — Spi): 

[cbso{x),M,%,iy)] = ^5^^\x - y) (62) 

We can now contemplate phases in which fluctuations of 
the spin density and the spin current are clamped. This 



is achieved by condensing vortices in (pso ■ Repeating the 
analysis of the previous section we see that the resulting 
paramagnet has (after refermionization) two species of 
fermions d„ at the original Fermi surface that are cou- 
pled to an emergent U{\) gauge field with opposite gauge 
charges: 



C — da \ dr — icao 



(V - iaay 
2m 



(63) 



This is an example of a spin-gapped "algebraic charge 
liquid" metal22i analogous to the holon metal that has 
been discussed in the context of doped Mott insulators. 
The coupling to the gauge field ultimate drives a pairing 
instability and the ground state will be superconducting. 
Nevertheless the spin gapped metallic ground state can 
be exposed if a magnetic field that suppresses supercon- 
ductivity is present. 



VII. HIGHER ANGULAR MOMENTA VORTEX 
CONDENSATES? 

The theoretical aproach described in this paper to 
the Mott transition immediately invites the question on 
whether more exotic phases may be accessed by con- 
densing vortices in higher angular momentum channels. 
Could we for instance introduce and condense vortices 
in the d^^-y^ channel of the charge boson, i.e. in 



f>c.2 + 0c,-2? The field 0c2 



_2 IS conjugate to a 



d,j.2_y2 distortion of the Fermi surface. Such a vortex 
condensate would have rather exotic properties: it will be 
a compressible metal which is nevertheless "incompress- 
ible" for such a d.j;2_y2 Fermi surface distortion. Further- 
more the single particle spectrum will be gapped almost 
everywhere on the Fermi surface as adding an electron at 
a generic fc-point necessarily couples in with the gapped 
dx2-y2 Fermi surface distortion mode. 

Despite the appeal of describing such an exotic metal 
very simply it is not clear to us that it is legitimate to 
ever introduce such higher angular momentum vortices 
and study their condensation. The difficulty is as fol- 
lows. To talk about vortices in some field it is important 
that it be a periodic variable. This in turn means that the 
physical Hilbert space is such that the conjugate variable 
is quantized to be an integer multiple of some basic unit. 
For the I = charge or spin boson the conjugate vari- 
ables are the total charge and spin densities respectively. 
Their spatial integrals can of course only take quantized 
values. Consequently 0co,0so are periodic variables and 
it is legitimate to study vortex defects. For the higher-^ 
modes however, the conjugate variables are various shape 
distortions of the Fermi surface which in the microscopic 
Hilbert space are not quantized. Consequently (pd or (psi 
for / 7^ should not be regarded as a periodic variable, 
and there is no obvious sense in defining vortices in these 
fields. 



VIII. SUMMARY 

In this paper we described a new approach to thinking 
about the electronic Mott transition on a two dimensional 
lattice between a Fermi liquid metal and a quantum spin 
liquid Mott insulator. We started with the effective low 
energy theory of the Fermi liquid viewed as a theory 
of shape fluctuations of the Fermi surface. The Mott 
transition must involve suppression of shape fluctuations 
that correspond to density and current fluctuations. To 
implement this suppression, and based on insights from 
the theory of the Mott transition of bosons, we incorpo- 
rated vortices in the electron phase into the bosonized 
description of the low energy description of a Fermi liq- 
uid. Condensing these vortices gives us a Mott insulator. 
We showed that the low energy effective theory of this 
Mott insulator has a spinon Fermi surface coupled to a 
fluctuating U{1) gauge field. This is the same result as in 
the usual slave particle description of the Mott transition 
even though our approach is very different and did not 
explicitly invoke slave particles. We also showed how we 
could obtain a description of a spin gapped metal with a 
charge Fermi surface by condensing vortices in the phase 
of the variable conjugate to the spin antisymmetric shape 
fluctuations of the Fermi surface. Such a state too has a 
well-known description in terms of slave particles. 

When compared to slave particle theories our approach 
is more coarse-grained. We keep only degrees of freedom 
near the Fermi surface and see how to change their dy- 
namics so that the Fermi liquid is destroyed in favor of 
a Mott insulator or a spin gapped metal. Thus we keep 
the momentum space structure of the Fermi liquid while 
incorporating the minimum possible correlation effects 
needed to describe the Mott insulator. In contrast the 
slave particle approach is very much a real space con- 
struct, and momentum space structure is introduced by 
first doing a mean field calculation to obtain a fermi sur- 
face for slave fermions. By providing an alternate point 
of view we hope that our approach gives additional in- 
sight into the novel physics associated with the Mott 
transition and the quantum spin liquid state. In partic- 
ular, our methods may provide a framework for thinking 
about phenomena such as the "momentum-selective Mott 
transitions" ^^'^° suggested by various cluster extensions 
of DMFT, or to describe exotic phases beyond what is 
accessible within a slave particle framework. We leave 
exploration of such possibilities for the future. 

We thank Tarun Grover and Brian Swingle for use- 
ful discussions. TS was supported by NSF Grant DMR- 
6922955. TS thanks the Aspen Center for Physics, where 
part of this manuscript was completed. 



Appendix A: Some calculations from the bosonized 
theory 

For clarity we will demonstrate the calculation of the 
commutator [(f)i{x),(l)i{y)] in the non-interacting case 



f{i,i') = 0. It should be clear however that the com- 
mutator is independent of f{i,i') since it is determined 
exclusively by the time-derivative terms in the action, 
while f{i,i') only appears with spatial derivatives of the 
fields. We compute 



= lim(0(a;, e), ^(y, 0) - ^{x, 0), 0(y, e)) 

c— >0 



(Al) 



— lim 



gifc-(x-y) gia;£ 



fcll,fc I 



AOkp fell (iw + VFk\\) 



-— hm - arctan — — (A2) 



. sgn(a:|| -y\i) D{xj_ - y±) 
2 AOkf ' 



(A3) 



where we employ a standard point-splitting procedure 
and the smeared delta- function D{x) is defined by the 
last line. Using Eq. ([T5|) we then obtain 



t n\ r^ An I nM ^ n/ \^ sgn(x||-y|,) 
i(a;,0),A^FA6ipi(j/,0)] = ^'^(xx ~yi_)d\\ ^ ^ 



27r 



(A4) 
(5(x|| -y||)D(a;x-y±). (A5) 



Next we are interested in the fermion correlator 

(V'(a;,r)V'^(0,0)) ~ exp47r2(0(a;, r)0(O, 0) - 0^^. (A6) 

In Ref. Il6| it has been shown that, in the long- wavelength 
limit, the fermion correlator is unaffected by a non- 
singular /(i,i'), so again we set it to zero for simplicity. 
It is then straightforward to compute 

(0(a;,T>O)(/.(O,O)-(/.2) (A7) 

-,icjO+ 



1 



Ax-h —\ru 



1 



2 



dk 



fc|| {ijj — \vpk\\\ 



gi2;||fc||g- 



(271")- jfe||>o 
In {LkpAe/TT) 

In (l — A|| {\x — w_ft)) 



for 



for 



xi_ > 



x±_ -C 



1 



kpAe 

1 
kpAe 



(A8) 
(A9) 



where L is the size of the system and Aj^ the momen- 
tum cut-off normal to the Fermi surface. We have now 
determined the fermion correlator 

{'4j{x,t)iIi\Q,Q)) ~exp47r2(0(a;,T)(/)(O,O) - </>2) (AlO) 

1 







^\' 



A|, — (ia; — vft) 



for 



for 



a;_L > 



x±_ « 



kpAe 

1 
kpAe' 



(All) 
(A12) 



where we have taken the thermodynamic limit L — )■ oo. 
Via the same point-splitting procedure as above, it is 
easily verified that the fermion operators satisfy anti- 
commutation relations. 
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Appendix B: Minimal coupling and vortices in 
bosonized action 



In the main text we point out that a bit of care is re- 
quired when couphng a gauge field to or isolating vortices 
from the bosonized theory. To illustrate the issue and its 
resolution, here we briefly discuss the 1-dimensional case, 
which may be more familiar to the reader. 

The bosonized action of one dimensional spinless 
fermions can be written in terms of the phases (J)r/l of 
the right/left moving fermions as 



S = Sr+Sl 



(Bl) 



Sr/L = T- I Tidr<pR/Ldx(f>R/L + {d.x<t>R/Lf- (B2) 

To couple the bosons to an electromagnetic field one may 
be tempted to follow the minimal substitution prescrip- 
tion 



5^0fl/L -^ 3p0fl/L + ^p 



(B3) 



This is however incorrect as we can easily see, e.g. by 
calculating the charge-density operator 



Pnaivc = i-TT- {SR{d^4>R + A^j) + Shid^^i'L + A^j)) 



dx{(t)B. -</'L)/47r, 



(B4) 
(B5) 



while we know that the correct charge and current oper- 
ators are given by 



/3 = dxi(t)R " (/>L)/27r 
J^dx{(bR + (j)L)/27r. 



(B6) 

(B7) 



This is not surprising. In fact, it is well known t hat the 
separation between right- and left-movers (Eq. (JB1|) ') is 
incomplete in the presence of gauge-fields due to the chi- 
ral anomaly^. Knowing the physical charge- and current 
operators we can readily write the correct action as 



S =SR{d^(l>R) + SL{d^(l>L) + 



/ -ipAr 

J X.T 



+ JA^ + ^ (B8) 

TT 



=SR{d^(pR + A^,) + SL{dfj,(t)L + A^ 



+ i{dr4>R + A.r)id^<l)L + A^)-{R^L). (B9) 



Note that the last term, which mixes left- and right- 
movers, vanishes upon integration by parts for A^ — 
and can then have an arbitrary coefficient. However for 
Af^ y^ the coefficient is uniquely given by Eq. (|B9I) . It 
is only for this choice of coefficient that minimal substi- 
tution yields the correct action. 
Equivalenty we may rewrite the action in terms of 



ip — -t) — 



(BIO) 



as 



5=2^ / -2idr^dx^ + {dxv^f + {dx^)^ (Bll) 



2tt 



-2i(a,</. + Ar)dx^ + {dxip + Axf + {dx^f. 

(B12) 



As above, for Ay, = 0, we are free to add 
/ {drfdx'& — drfdx'd) to Eq. (|B11[) with a arbitray co- 
efficient, but minimal substitution only yields the correct 
answer (jB12[) if the coefficient is zero. 
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